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Abstract 

We use the Nambu-Jona-Lasinio model as an effective quark theory 
to investigate the medium modifications of the nucleon electromag- 
netic form factors. By using the equation of state of nuclear matter 
derived in this model, we discuss the results based on the naive quark 
- scalar diquark picture, the effects of finite diquark size, and the me- 
son cloud around the constituent quarks. We apply this description to 
the longitudinal response function for quasielastic electron scattering. 
RPA correlations, based on the nucleon- nucleon interaction derived in 
the same model, are also taken into account in the calculation of the 
response function. 
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1 Introduction 



The structure of the nucleon and its modifications in the nuclear medium is a 
very active field of experimental and theoretical research. The basic quanti- 
ties, which reflect the charge and current distributions in the nucleon, are the 
electromagnetic form factors pQ, which are currently investigated in elastic 
electron- nucleon scattering experiments from intermediate to very high en- 
ergies The knowledge of the nucleon form factors is also inevitable to un- 
derstand the electromagnetic structure of nuclei. Electron-nucleus scattering 
experiments under quasielastic kinematic conditions, like the measurement 
of inclusive response functions in the intermediate energy region E] and 
recent measurements of polarization transfer in semi-exclusive knock-out pro- 
cesses j3] , are ideal places to study the form factors of a nucleon bound in the 
nuclear medium. Because the structure of the quark core and the surround- 
ing meson cloud may be different for a bound nucleon and a free nucleon, one 
expects medium modifications of the nucleon form factors [6 , and the explo- 
ration of these effects is an important subject at current electron accelerator 
facilities [7j. 

On the theoretical side, effective quark theories are the ideal tools to 
describe the electromagnetic form factors of the nucleon. Much progress for 
the case of the free nucleon has been made in Faddeev type descriptions based 
on the Schwinger- Dyson method jH]. An important point which still has to 
be implemented in these calculations is the role of the pion cloud around the 
nucleon, and the recently developed method of chiral extrapolations of lattice 
results jHj provides important hints. On the other hand, the calculation 
of form factors at finite nucleon density requires also a description of the 
equation of state of the many-nucleon system, and here progress has been 
made by using the Nambu-Jona-Lasinio (NJL) model JU] as an effective 
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quark theory: Recent works have shown how to account for the saturation 
properties of nuclear matter in this model and when combined with the 
quark-diquark description of the single nucleon^] this provides a successful 
description of both nucleon and nuclear structure functions for deep inelastic 
scattering [131 E] 1 . 

The purpose of this paper is to discuss the results for the nucleon form 
factors obtained in the simple quark-scalar diquark description of the nucleon 
at finite density in the NJL model. We have to note from the beginning that 
this can only be a first step toward a realistic description, because it is known 
that axial vector diquarks are important for spin- dependent quantities |17| 
IT4] . and the pion cloud is important for magnetic moments and the size of 
the nucleon CHI- While the axial diquarks could be included in a further 
step like it was done for the structure functions [14j, a reliable description 
of pion cloud effects makes it necessary to go beyond the standard ladder 
approximation scheme. However, like the simple quark-scalar diquark model 
of Ref . ^3] served as a basis for the more elaborate description of structure 
functions it will also be the basis of a more realistic description of form 
factors including axial diquarks and the pion cloud. To provide this basis is 
the main intention of the present paper. 

In Sect. 2 we will briefly review the model for the nucleon and the nuclear 
matter equation of state. Sect. 3 is devoted to the nucleon form factors at 
finite density, and in Sect. 4 we discuss the numerical results. As an applica- 
tion, we discuss the response function for quasielastic electron scattering in 
Sect. 5. For this purpose we will also elucidate the nucleon- nucleon interac- 
tion in our model in order to include the correlations within the relativistic 
RPA. A summary will be presented in Sect. 6. 

1 Recently the model has been extended to describe the equation of state at high 
densities [T51 IT?)] . 
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2 The model 



In this work, we use the NJL model as an effective quark theory to describe 
the nucleon as a quark-diquark bound state, and nuclear matter (NM) in the 
mean field approximation. The details are explained in Refs. fUEI]) and 
here we will only briefly summarize those points which will be needed for our 
calculations. 

The NJL model is characterized by a chirally symmetric 4-fermi interac- 
tion between the quarks [T^j. Any such interaction can be Fierz symmetrized 
and decomposed into various qq channels [201 • Writing out explicitly only 
those channels which are relevant for our present discussion, we have 

C = ^(^-m)^ + G 7r ((W) 2 - (V(7 5 r)^) 2 ) - tyftf + • • • (2-1) 

where m is the current quark mass. In a mean field description of the isospin 
symmetric nuclear matter ground state \p), the Lagrangian can be expressed 
as 

C = W-M-y)*-&£& + V£ + C,, (2.2, 

where M = m — 2G 7r (p\ipilj\p) and = 2G w (p\ip'y ll ifj\p), and d is the 

normal ordered interaction Lagrangian. The effect of the mean scalar field is 

thus included in the density-dependent constituent quark mass M, and the 

effect of the mean vector field is to shift the quark momentum according to 

P^ = Pq + ^ A '> where Pq is the kinetic momentum. The propagator of the 

constituent quark therefore has the following dependence on the mean vector 

field 2 : S{k) = S{k Q ). 

2 In this section, Green functions in the presence of the mean vector field are denoted 
by a tilde, and those without the vector field have no tilde. In the loop integrals for 
the electromagnetic form factors in Sect. 3, however, it is always possible to eliminate the 
vector field by a shift of the integration variable, and therefore the tilde-Green functions 
do not appear in later sections. 
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One can use a further Fierz transformation to decompose Cj into a sum 
of qq channel interaction terms [201 • For our purposes we need only the 
interaction in the scalar diquark ( J n = + , T = 0, color 3) channel: 

Cj, s = G s ( 75 C) r 2 (3 A ^ T ) (> T (C7- 1 T5 ) r 2 (3 A i) , (2.3) 



where (3 A = y3/2 X A (A = 2,5,7) are the color 3 matrices and C = ijzjo- 
The coupling constant G s will be determined so as to reproduce the free 
nucleon mass. 

The reduced t-matrix in the scalar diquark channel is given by jTTj 

AiG 

fJq) = -^4 = T s (q B ) (2.4) 

with the scalar qq bubble graph 

= 6i / 7^ tr ^ [75^(^)75^ (-(<? - fc))] = IT s (g D ) . (2.5) 

Here = q^ — 2 V 1 is the kinetic momentum of the diquark. 

The relativistic Faddeev equation in the NJL model can been solved nu- 
merically for the free nucleonj20], but here we restrict ourselves to the static 
approximation [21], where the momentum dependence of the quark exchange 
kernel is neglected. The solution for the quark-diquark t-matrix then takes 
the simple analytic form 

with the quark-diquark bubble graph given by 

n*(p) = - / -J^ s(k) u v -k) = u N ( PN ) , (2.7) 

where p^ = p M — 3V 1 is the kinetic momentum of the nucleon. The nucleon 
mass Mm is defined as the pole of ()2.6|) at ^ N = Mn, and the positive energy 



spectrum has the form p = e p = E Np + 3V with E Np = \J Mjj + p 2 . The 
nucleon vertex function in the non-covariant normalization is defined by the 
pole behavior of the quark-diquark t-matrix: 

r N ( P )T N {p) 

T N {P) -> as p -> e p . (2.8) 

Po ~ e P 

From this definition and Eq. ()2.6|) one obtains 



^n(p) = \ ~Z N — ^ u N (p N ) , (2.9) 

V &Np 

where un is a free Dirac spinor for mass Mm normalized as un^n — 1- The 
normalization factor is easily obtained from this definition and will be 
given in Eq. (|2.14jl below. We note that with this normalization the vertex 
function satisfies the relation 

In the numerical calculations of this paper, we will approximate the quantity 
r s by a "contact+pole" form: 

T s (q) -> AiG s - 19 8 = AiG s - ig s A Fs {q). (2.11) 
q z — M£ 

Here A^ s (g) is the Feynman propagator for a scalar particle of mass M s , 
which is defined as the pole of r s of Eq. (j2.4j) . The residue at the pole (g s ) 
will be given in Eq. ()2.13j) below. 

In the calculation of the nucleon form factors, we will also consider the 
effects of the pion cloud around the constituent quarks. In this case, the 
propagator S(p) of the quark involves an additional self energy correction 
from the pion cloud (Sq). Here we will use a simple pole approximation for 
S(p): 

S(p) = Z Q S F (p), Z Q l = l-(^) , (2.12) 



dp" 



=A1 
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where S F is the Feynman propagator of a constituent quark with mass M. 
In this approximation, the pion effects can be renormalized by i/j — > \J Zqip 
and a redefinition of the four fermi coupling constants G a — > Ga/Zg, see 
Ref.( 12J). For the calculation of the form factors, however, we will keep 
the factor Zq explicitly for clarity 3 . Introducing (J2.12)) and (J2.ll)) into the 
expressions (12 .5j) and (|2.7|) for the bubble graphs shows that the diquark and 
nucleon normalization factors can be written as 



i t aa,( q ) \ i yi I dn.( q ) . _ . , 



(2.13) 



^n 1 — tt, — —Z Q g s — — = ZqggZj^ 1 , 

V )^=M N \ W )i>=M N 

(2.14) 

where fl s and IT at are defined in terms of the pole parts only: 

ft.(ff) = QiJ^- 4 tY D [ l5 S F (k) l5 S F (-(q-k))\ (2.15) 
fl N (p) = ij^S F (k)A Fs (p-k) (2.16) 

The equation of state of NM in the NJL model can be derived in a formal 
way |m from the quark Lagrangian ()2.2|) by using hadronization techniques, 
but in the mean field approximation the resulting energy density of isospin 
symmetric NM has the simple form JT] 

£ ^-^ + 4 /|^-|pI)^ < 2 - 17 ' 

where the vacuum contribution (quark loop) is 

, d 4 A; fc 2 -M 2 + *e (M-m) 2 (M - m) 2 
gv = 12 V (27F ln A: 2 — Mq + ie + ~^ 4Gr~ " (2 ' 18) 



3 We also note that such a renormalization procedure is no longer possible when one 
considers the pion cloud effects around the nucleon, which goes beyond the simple ladder 
approximation on the quark level. 
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Here M the constituent quark mass for zero nucleon density. The condition 
dS/dVo = leads to V Q = 6G w p, and we can eliminate the vector field in 
(I2.17|) in favor of the baryon density. The resulting expression has then to be 
minimized with respect to the constituent quark mass M for fixed density. 
For zero density this condition becomes identical to the familiar gap equation 
of the NJL model^U], and for finite density the nonlinear M-dependence of 
the nucleon mass Mn is essential to obtain saturation of the binding energy 
per nucleon [TT] . 

In order to fully define the model one has to specify a cut-off procedure. 
In the calculations in this paper we will use the proper time regularization 
scheme |22[ lllj. where one evaluates loop integrals over a product of propa- 
gators by introducing Feynman parameters, performing a Wick rotation and 
replacing the denominator (A) of the loop integral according to 



where Air and Auv are the infrared and ultraviolet cut-offs, respectively. 
The infrared cut-off plays the important role of eliminating the unphysical 
thresholds for the decay of the nucleon and mesons into quarks [22] , thereby 
taking into consideration a particular aspect of confinement physics. 

3 Nucleon electromagnetic form factors 

The electromagnetic current of the nucleon in the quark-diquark model is 
represented by the Feynman diagrams of Fig.l and given by 



Here Aq and A^, are the electromagnetic vertices of the quark and the 
scalar diquark, both depending on the final and initial particle momenta. It 




(n > 1) 



(2.19) 




+i {r s (p' - k)A» D T s (p - k)) S{k)\ T N (p) . 



(3.1) 
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Figure V. Feynman graphs for the evaluation of the nucleon electromagnetic 
current in the quark-diquark model. The double line represents the diquark 
t-matrix, the solid line the constituent quark propagator, and the black areas 
are the vertex functions and electromagnetic vertices. 

is an easy task to use the Ward-Takahashi identities [23! f° r the Q uar k and 
diquark vertices 

qtl (s(£'W Q S(£)) = -Q Q (S(£')-S(£)) (3.2) 
?,(r s (f)A^r s (f)) = iQ d (t s (£')-t s (£)) (3.3) 

to show the current and charge conservation for the nucleon: 

qMq) = QNT N (p')(n N (p')-n N (p))T N ( P ) = o (3.4) 

f(0) = Qnj^, (3.5) 

where the electric charge of the nucleon Q N = Qq + Q D is the sum of the 
quark and diquark electric charges. 

The electromagnetic vertices in ()3.1|) describe the finite extension of the 
constituent quarks and the diquark. In general, they should be calculated 
off-shell consistently with the propagators r s and S, using Feynman diagrams 
or some ansatz which satisfies the Ward identities[H]. In the present work, 
however, our principal aim is to investigate the medium effects in a simple 
model calculation. For this purpose, we limit the complications caused by 
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the quark and diquark sizes to a minimum, and use an on-shell (or pole) 
approximation for the quark and diquark currents appearing in Eq. (|3.ip : 

(S(e')A^S(l)) — > Z Q (s F {l')k» Q S F {l)) (3.6) 
(r s (f)A;r s (f)) — - S ,(A ft (OA£A Fl W), (3.7) 

where the on-shell (o.s.) vertices are denoted by a hat and given by the pole 
residues of the full quantities by 

A£ ee Z Q (A» Q ) os =YF 1Q (q 2 ) + ^F 2Q tf) (3.8) 

a£ ss ^(Aa.,^'+^r^(g 2 ). (3.9) 

Here we introduced the quark and diquark form factors which satisfy -Fiq(O) = 
Q Q , F d (0) = Q D . 

To understand (|3.6|) and (|3.7|) . we note that in general the on-shell approx- 
imation for a vertex function A M (£',£) can be formulated only if it appears 
between pole parts of Green functions, because only in this case one can ap- 
proximate it by its value for on-shell momenta £', I 4 . This is why in Eq. fl3.6j) 
and (J3.7)) we have replaced the propagators left and right to the vertex func- 
tions by their pole parts (see (j2.11|) and (|2.12j) ). which is also essential in 
order to have charge conservation with the vertices ()3.8|) and ()3.9|) . 

We now can write down the form of the nucleon current (|3.1|) which will 
be used in the further calculations: 

f N (q) = (0£ + 0% + 0£) u N (p) . (3.10) 

Here the first and second terms denote the contributions of the contact term 

and the pole term of the diquark t-matrix to the quark diagram (first diagram 

4 This approximation is the basis of the standard convolution formalism to calculate 
quark light-cone momentum distributions and structure functions. 
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of Fig.l), and the third term is the contribution from the diquark diagram: 

O'c = ~^Z N J -^S F (p> - *) (^Q(^) + W^^g(g 2 )) S F (P ~ k) 

(3.11) 

0% = iZ N J j^S F (j/ - k) (i»F 1Q (q 2 ) + ^iW)) S F (P ~ k)A Fs (k) 

(3.12) 

0£ = iZ N F D (q 2 ) J -^^Fs(p' -k)(p + p'- 2kf A Fs (p - k) S F (k) . (3.13) 

(For the contact term, we replaced G s — ► G s /Zq, so that G s in (|3.11|) is 
the renormalized coupling in the sense explained in Sect. 2.) By using the 
elementary Ward-Takahashi identities = Sp 1 {(!.') — Sp {£) and £' 2 — £ 2 = 
A Fs (£') — A Fs (£) and the fact that on the nucleon mass shell H-n{p) = H-n(p'), 
it is easy to check that the 3 parts in (|3.1U|) satisfy current conservation 
separately 5 . Similarly, charge conservation can be checked by using the el- 
ementary Ward identities Y = dSp l {£)/d£^ and W = dA F \{£)/d£^, as 
well as F\q(0) = Qq, -Fd(O) = Qd- It has to be noted, however, the general 
Ward-Takahashi identity for an off-shell nucleon (the first equality in Eq. ()3.4|) 
without the nucleon spinor) is not valid in this approximation scheme. 

It is not very difficult to evaluate the three loop integrals in (|3.11|) - (|3.13|) . 
and in Appendix A the results are given in terms of the Dirac-Pauli form 
factors Fin and F 2 n, which are defined by 

IM N M N _ , r 2 i^q v . 2 J , s 

j n = \ TT- ~B — u n(p ) l^F 1N {q ) + ——F 2N {q ) u N {p) . 

(3.14) 

In the following we will discuss various steps in the calculation of the 
nucleon form factors. 

5 These formal manipulations involve shifts of the integration variables. In the actual 
calculations based on our regularization scheme, however, we always checked that current 
and charge conservation are satisfied exactly. (Therefore, the explicit expressions given in 
the Appendices all satisfy the charge conservation.) 
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3.1 Naive quark-diquark model 

The simplest approximation consists in assuming point couplings of the 

1 r 3 

quarks and diquarks to the photon, i.e., replacing Fin — > Qn = — + — , F 2 o — > 0, 

""62" 

Fd — > = - in Eqs. (j3.11j) - (j3.13j) . This approximation will be called the 
"naive quark-diquark model", and the detailed expressions can be found in 
Appendix A. 

3.2 Effects of finite diquark size 




Figure 2: Graphical representation of the diquark form factor. The white 
circles denote the diquark vertex functions, and the black area is the quark 
electromagnetic vertex. (There is a second diagram where the photon couples 
to the other quark, but also an overall symmetry factor |.) 



Here we consider the effect of the diquark form factor Fjj, which has been 
defined in ()3.9|) . The vertex A^, is shown graphically in Fig. 2, where the 
quark-diquark vertex functions are those appearing in the Lagrangian (j2.3|) . 
We obtain 

A» D = ig s J^{Tr{ l5 S(pi + k)A» Q S(p+k) l5 S(k)}} , (3.15) 

where the trace refers to color, isospin and Dirac indices. Using the on- 
shell approximation ([3.6)1 . the definition of quark form factors ([3.8)1 . and also 
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(gH2D and ([HTgD , we obtain 



AM 



6i g s 



(2tt)< 



+ F$(Q 2 )Tr D 



{F^(Q 2 )Tr D [ l5 S F (p' + k)^S F (p + k) l5 S F (k)} 

(3.16) 



lbS F (p' + fc) nn / S F (p + k)j 5 S F (k) 



2M 



Here Fyj and F^q are the isoscalar parts of the quark form factors Fig and 
F 2 q. The resulting diquark form factor Fjj is given in Appendix A. 

3.3 Effects of meson cloud around constituent quarks 

Here we consider the quark form factors arising from the pion cloud and 
vector mesons. For the pion cloud, we obtain from the definition (|3.8J) and 
the Feynman diagrams of Fig. 3, 

d A k 



A£ = Z Q rQ Q + Z Q 



(2vr) 4 

+ i( T7r (p'-k)A%T n (p-k)) l5 S(k) l5 } 



7s Ti (S(p' - k)A^ S(p - kfj TasT^k) 

(3.17) 



i \ 




Figure 3: Pion cloud contributions to the quark electromagnetic vertex. Zq 
is the quark wave function renormalization factor, and the dashed line rep- 
resents the t-matrix in the pion channel. (This figure actually refers to the 
expression ([3.25)1 . where the factor Zq remains only for the "bare" term.) 



Here t w is the reduced qq t-matrix is the pion channel, which can be 
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approximated in the same way as the diquark t-matrix (J2.11j) : 

T«(k) = i - ~^ {k2) — -XG« + = -2iG. + ig^{k) . 

(3.18) 

The bubble graph in the pion channel is 11^ = IT S (see 1)2. 5|) ). and the residue 
g w will be given in Eq. (|3.24j) below. 

The expression (13.17)) is formally similar to the nucleon current (|3.1|) . The 
presence of the first term simply expresses the fact that in the NJL model 
the (bare) quarks are present from the beginning, and the factor Zq gives 
the probability of having a constituent quark without its pion cloud. It has 
been expressed in Eq. (|2.12j) in terms of the self energy 

£g(p) = "3 / (lsS(p - % 5 ) r K (k) . (3.19) 

The quark electromagnetic vertex Aq Q in ()3.17|) will be approximated by 
its point form after processing the renormalization factor Zq, and the pion 
electromagnetic vertex = TiA^^Tj is similar to the diquark vertex of Fig.2, 
but with point quark-photon couplings, as will be specified below. 

We now follow the same steps as for the calculation of the nucleon current: 
We use the on-shell approximation for the quark and pion vertices 

(S(£')A» Q0 S(£)) — > Z Q (S F (£')A% S F (£)) (3.20) 
(r w (£')A^(£)) — -g w (A Fw {£')A^A Fw (£)) , (3.21) 

where the on-shell vertices are defined by 

A£ = 2g(A5o) M =YQ Q (3.22) 
= g*(K#) M = (£' + £T He ij3 ) F^q 2 ) . (3.23) 
Using S(k) = ZqSpik) in the expression for U n , we get 
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with the renormalized bubble graph II 



-7T 



n a , see (jZTH|) . Then Eq. (|3~T7T 



becomes 




Zq^Qq - ig,\{\ - r 3 ) / (S F (p' - k)>fS F (p - fc)) 7 5 A^(fc) 

22r 3 ^^(g 2 ) | A^(A Fw (p , -A;)(p / + p-2A ; ) M A^(p-A ; ))7 5 5 F (fc)75. 

(3.25) 



Here we note that the contribution of the contact term (2iG^) to the quark 
diagram has been dropped in order to avoid double counting: Because we 
always assume that our interaction Lagrangians are Fierz symmetric, it is 
easy to see that this contribution can be incorporated into the vector meson 
channel, which will be separately considered below. By using S — Zq S f in 
the self energy (j3.19|) and in the expression for Zq of Eq. (j2.12j) . we see that 



The further evaluation of the loop integral ()3.25|) is left to Appendix B, where 
the contributions to the quark form factors F\q and F2Q are given. The pion 
electromagnetic vertex is evaluated from the definition ()3.23|) and a Feynman 
diagram similar to Fig. 2 for an external pion, but with a point quark-photon 
coupling: 




(3.26) 



where the renormalized quark self energy is given by 




(3.27) 




I 



d 4 k 



Tr D [ 75 S F (p' + k)<y"S F (p + k) 75 S F (k)} . (3.28) 



The explicit form of F^{q 2 ) is given in Appendix B. 
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Finally, we consider the corrections of the quark electromagnetic vertex 
arising from vector mesons, similar to vector meson dominance (VMD) mod- 
els. If our original interaction Lagrangian contains terms of the form 

Ci, v = -G u (?7^) 2 - G p (^r^) 2 , (3.29) 



then the point-like quark-photon vertices in the diagrams of Fig. 3 and the 
pion vertex are replaced by the VMD vertex shown in Fig.4. 




Figure 4: The corrections to the quark electromagnetic vertex arising from 
vector meson dominance. Here the 4-fermi interactions refer to the vector 
channel (see Eq. (|3.29|) ). i.e., terms with tensor coupling to the external quark 
are not included. 



Because of the transverse structure of the bubble graphs in the vector 
channel, this leads to the following renormalization: 



+ 



1 

6 

2 



r 



r 



2Gji v {q 2 ) 

l + 2G ( JL / (g 2 ; 
2G p Il v (q 2 ) 

l + 2G p fl v (q 2 



r 



r 



q 2 ) 

I 2 j 



(3.30) 
(3.31) 



where the form of Ily is given in Appendix B. Because our quark-photon ver- 
tex in ()3.8p is defined for on-shell quarks, the terms cx — — do not contribute. 
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Therefore the isoscalar (or isovector) parts in the quark electromagnetic ver- 
tex should be multiplied by a form factor 6 F^q 2 ) (or F p (q 2 )), where 

F «tf) = - - or \ . 2 . {cl = u,p). (3.32) 
1 + 2G a Tl v (q 2 ) 

4 Results for the nucleon form factors 

In this section we will show our results for the nucleon form factors. First 
we discuss our model parameters and particle masses. The parameters are 
the same as in Refs. |ll| ITS] : The IR cut-off Air is fixed as 0.2 GeV, the 
constituent quark mass at zero nucleon density is M = 0.4 GeV, and G n , 
A uv are determined so as to reproduce M w = 0.14 GeV and f w = 93 MeV. 
This gives G n = 19.60 GeV~ 2 and Auv — 0.6385 GeV. The coupling constant 
G s is determined so as to reproduce M^o = 0.94 GeV, which gives the ratio 
G s /G n = 0.508. The coupling constant G^ is determined so that the curve 
for the NM binding energy per nucleon (Eb/A) as a function of the density 
passes through the empirical saturation point 7 (p,Eb/A) = (0.16 fm -3 , 15 
MeV), which gives the ratio G^/G^ = 0.37. Finally, for the VMD form 
factors ()3.32|) we also need the coupling constant G p , which is determined by 
reproducing the empirical symmetry energy coefficient = 35 MeV at the 
density p = 0.16 fm -3 , which gives G p /G n = 0.091 

In Table 1 we list the effective quark, diquark, nucleon and pion masses 
for the densities p = 0, 0.08, 0.16, and 0.24 fm -3 . Concerning the pion mass 
in the medium, we use a general result based on chiral symmetry [23] , which 
for the NJL model implies that the product M 2 -M is a constant independent 

6 Because the scalar diquark has isospin zero, this eventually also holds for the nucleon 
electromagnetic vertex. 

7 We recall from R.ef. jllj that, in this simple NJL model, we cannot adjust both the 
empirical binding energy and saturation density at the same time. Therefore, although 
the binding energy curve passes through the empirical saturation point, its minimum is at 
a different point, (p,E B /A) = (0.22 far 3 , 17.3 MeV). 
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of density, see Eq.(2.58) of Ref.[TT]- We therefore use 8 Ml 
Also listed in Table 1 are the values of g s , g w , Z^ and Zq. 



M* ■ M /M. 





p = 


p = 0.08 fm- 3 


p = 0.16 fm 3 


p = 0.24 fm~ 3 


M 


0.4 


0.353 


0.308 


0.268 


M s 


0.576 


0.493 


0.413 


0.342 


M N 


0.94 


0.818 


0.707 


0.619 


M n 


0.14 


0.149 


0.159 


0.171 


9s 


18.26 


17.17 


16.20 


15.40 


£h 


17.81 


14.29 


11.60 


9.70 


Z N 


44.32 


47.62 


50.38 


51.97 


Zq 


0.808 


0.832 


0.852 


0.869 



Table 1: Effective masses for the quark, the diquark, the nucleon and the 
pion (all in GeV), and pole residues g s , g T , Z N and Zq for four values of the 
density. 

We will discuss our results for the nucleon form factor in terms of the 
Dirac-Pauli form factors defined by Eq. (|3.14|) . In the discussion of medium 
effects, in particular the effects of the reduced nucleon mass (M N < M N0 ), 
we will also refer to an equivalent parametrization in terms of the "orbital 
form factor" Gl and the "spin form factor" Gs'- 

(p' + pY 
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2M 
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-G L (q 2 ) + 
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Gs(q 2 ) 



u N (p) 



(4.1) 



The relations to the Dirac-Pauli form factors are 
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2N, 



(4.2) 



8 More precisely, this pion mass is defined at zero momentum, and includes nucleonic 
(Z-graph and contact) terms, which are important in order to guarantee the Goldstone 
nature of the pion in the medium. We note that these nucleonic contributions to the 
scalar diquark (or sigma) mass at normal densities are numerically small compared to the 
qq (or qq) polarizations, although they become important for small M and guarantee the 
stability of the system w.r.t. variations in M, see Ref.(^J) for details. 
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We introduce these form factors here, because F\ and F 2 do not directly 
reflect the enhancement of the nucleon orbital current (j N orb ) arising from 
the reduced nucleon mass (enhanced nuclear magneton) 9 . Moreover, the 
parametrization (|4.1j) for the space part of the current has more connection 
to the traditional calculations of nuclear magnetic properties j2H] 5 because 
the values of these form factors at q = reduce to the orbital and spin 
g-factors: g L = G L (0), g s = G s (0). 

Here we would like to point out that these different ways to discuss the 
medium modifications of the nucleon current remind us that the form factors 
of a nucleon in the medium are not directly observable quantities: Ultimately 
the current has to be used in a nuclear structure calculation of observable 
cross sections. Our current reflects only those effects which are not taken 
into account in nuclear structure calculations, i.e., the effects of the nuclear 
mean fields on the internal motion of quarks in the nucleon. Other effects, 
which explicitly depend on the density and have their origin in the Pauli 
principle on the level of nucleons, must be considered in the nuclear part 
of the calculation 10 . As an example of such a calculation for the case of 
nuclear matter, we will consider the response function for quasielastic electron 
scattering in Sect. 5. 

For the zero density (single nucleon) case, it is possible to define a Breit 
frame where go — 0, and in this frame the nucleon current can be expressed 
in terms of the familiar electric and magnetic form factors 

G E (q 2 ) = Fi(g 2 ) + ^F 2 (q 2 ) (4.3) 

G M {q 2 ) = F^q 2 ) + F 2 (q 2 ) . (4.4) 

9 Note that the appearance of the medium modified nucleon mass (-Mjv) in the Pauli 
term of Ij3.14|) is a mere definition of Fi . 

10 This is also evident from the fact that the full current of a nucleon in the medium, 
including the explicitly density dependent parts, cannot be parametrized in the Lorentz 
invariant forms i|3.14[l or l|4.1|l . 
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We will compare our calculated form factors for zero density with the empir- 
ical dipole form factors, defined by 

Gep — Gupj ^p — Gmu/ l-ln — ~, 

(l - g 2 /0.71 GeV 2 ) 
F ln = 0. (4.5) 

For a nucleon moving in the medium, however, one cannot define a Breit 
frame, and consequently the combinations ()4.3j) and (|4.4|) do not enter nat- 
urally in the expressions for nuclear observables, like response functions or 
elastic form factors. For the finite density case we will therefore discuss our 
results in terms of the form factors F\ and F 2 , or Gl and G$- 

The results for the form factors at zero density (free nucleon case) are 
shown in Figs. 5-8. There we plot (i) the results of the naive quark-diquark 
model (see Sect. 3.1) without (dotted lines) and with (dashed lines) the con- 
tact term contribution to the quark diagram of Fig.l, (ii) the results obtained 
by including in addition the effects of the diquark form factor (dash-dotted 
lines), and (iii) the total result including also the pion and VMD effects. 

Fig. 5 shows that in the naive quark-diquark model the electric size of the 
proton is too small and the form factors F lp and Ge p fall off too slowly. The 
situation improves when the diquark form factor is included, and also the 
pion cloud gives some positive contribution to the electric size of the proton 
11 . The total result for F\ p still lies above the empirical dipole form factor, 
but we can expect that the further inclusion of axial vector diquarks and 
pion cloud effects around the nucleon will lead to a satisfactory description. 

Fig. 6 shows that the proton magnetic moment in the naive quark-diquark 
picture is too small, which is expected and also well known [TTl l24|. The finite 

n We obtain < r% > p = 0.421 fm 2 + 0.062 fm 2 = 0.483 fm 2 , where the two terms come 
from the Dirac (F±) part and the anomalous (F2) part. The fact that this is too small 
compared to the experimental value of 0.74 fm 2 is partially because the magnetic moment 
is too small, but also because the slope of Fi is too small. 
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Figure 5: The form factor F\ p for zero density. The dotted line is the result 
of the naive quark-diquark model without the contact term, to which the 
other contributions are successively added as follows: Dashed line: including 
contact term; dash-dotted line: including diquark form factor; solid line: 
full results including pion and vector meson contributions. The dash-double 
dotted line labeled "dipole" shows the dipole parametrization. The insert 
shows the form factor Ge p , with the same meaning of the lines. 

size effects of the scalar diquark do not contribute much in this case. The 

inclusion of the contact term in the quark diagram of Fig.l improves the 

situation somewhat. By using a Fierz transformation, this term is actually 

seen to be equivalent to a vector meson contribution (like shown in Fig.4 for 

the quark), but with a tensor coupling (oc cr^ u q u ) to the nucleon. Also the 

pion cloud, which leads to anomalous magnetic moments of the constituent 

up and down quarks 12 , gives a positive contribution to the proton magnetic 

moment, but the total result is still too small. It is, however, known that the 

12 We obtain fx u = | + 0.061, fid = — | — 0.123 for the magnetic moments of u,d quarks 
in the free nucleon. 
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0.0 0.5 „ 1.0 „ 1.5 2.0 

Q [GeV 2 ] 

Figure 6: Same as Fig. 5 for the form factors Fi v (main part) and Gmp 
(insert) for zero density. 



axial vector diquark and the pion cloud around the nucleon give large contri- 
butions to the magnetic moment and the associated form factors, and Fig. 6 
shows how far one can go in the simple quark - scalar diquark description. 

The importance of the diquark form factor is also seen for the neutron 
form factors F\ n and Geu, which are shown in Fig. 7. The naive quark- 
diquark model gives an electric form factor which is too large in comparison 
to the experimental one (note that the experimental Geu is smaller than 
the "dipole form factor" shown in Fig. 7), and the diquark form factor, 
which suppresses the (positive) contribution from the second diagram of Fig. 1 
relative to the (negative) first one, is essential to obtain reasonable values. 
The electric size of the neutron is somewhat too small in magnitude, but 
this is because the absolute value of the magnetic moment is too small 13 . 
13 We obtain < r| > n = 0.003 fm 2 - 0.072 fm 2 = -0.069 fm 2 , where the two terms come 



22 



In this connection, it is interesting to observe that the result for and in 
particular its contribution to the electric radius, is very small, and therefore 
the electric size of the neutron is almost entirely due to the "Foldy term" 

0.20 I 1 1 r 

0.15 — 



0.10 — 



0.05 — 
0.00 £77 



-0.05 

0.0 0.5 „ 1.0 „ 1.5 2.0 

Q [GeV 2 ] 

Figure 7: Same as Fig. 5 for the form factors F\ n (main part) and Geu (in- 
sert) for zero density. The dipole parametrization of F\ n is zero by definition 
and therefore not indicated here. 




The situation for the form factors i<2 n and Gmu shown in Fig. 8 is similar 
to the case of the proton (Fig. 6), i.e., the contact term and the pion cloud 
around the constituent quarks give some improvements of the magnetic mo- 
ment, but the total result is still too small. This, and the fact that the form 
factors fall off too slowly, again points out the necessity to include the axial 

from the Dirac (F±) part and the anomalous (F2) part. The fact that this is too small in 
magnitude compared to the experimental value of — 0.12fm 2 is because the absolute value 
of the magnetic moment is too small. For discussions on the role of these two contributions 
to the neutron electric radius, see for example 26 . 



23 



vector diquark channel, 
o.o 
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Figure 8: Same as Fig. 5 for the form factors F 2n (main part) and Gmu 
(insert) for zero density. 

The medium modifications of the nucleon form factors are shown in Figs. 
9-12, where we plot the results for p = (dotted lines), p = 0.08 fm -3 (dashed 
lines), p = 0.16 fm -3 (solid lines), and p = 0.24 fm -3 (dash-dotted lines). 

The result for F lp of Fig. 9 indicates that the electric size of the proton in 
the medium is somewhat enhanced. The orbital form factor Gl p shown in the 
insert of Fig.9 demonstrates the enhancement of the orbital current (jfNorb) 
due to the reduced effective nucleon mass, see Eq. (J4.2|) . We have to remind, 
however, that the isoscalar part of this enhancement is in a sense spurious, 
because in an actual nuclear calculation it is canceled by the "backflow" 
effect [27., which in our language arises from Z-graphs, i.e., the Pauli blocking 
part of the NN excitation piece (see the detailed discussions in Ref . [281 12*9*] 
on the backflow in relativistic meson- nucleon theories). Namely, the proton 
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Figure 9: The form factor F\ p for the cases p = (dotted line), p = 0.08 fm -3 
(dashed line), p = 0.16 fm -3 (solid line), p = 0.24 fm -3 (dash-dotted line). 
The insert shows the proton orbital form factor Gl p defined in Eq. (|4.1j) with 
the same meaning of the lines. 

orbital g-factor, which is roughly M^ /M^ in a Hartree calculation, becomes 
approximately |(1 + M^o/M^) after the inclusion of the backflow, where the 
first term is the isoscalar and the second one the isovector piece 14 . 

Fig. 10 shows that the medium effects tend to decrease the "intrinsic" 
anomalous magnetic moment of the proton, but when combined with the en- 
hancement of the nuclear magneton, the spin g-factor is enhanced, as shown 
in the insert of Fig. 10. It is interesting that a very similar result has been 
obtained also in hadronic mo dels [TH]. Therefore, the quark effects considered 
here do not lead to a quenching of the spin g-factor, as would be desirable to 
explain the missing quenching of isovector nuclear magnetic moments |2*S] . 
14 It is well known that the pion effects further enhance the isovector piece |28l. 
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Figure 10: Same as Fig.9 for the form factor i*2 P (main part) and the proton 
spin form factor G$ p defined in Eq. (|4.1|) (insert). 



but rather to an enhancement 15 . The figure also shows that the magnetic 

size of the proton becomes somewhat larger in the medium. 

The results for the neutron form factor F\ n in Fig. 11 show that the 

effect of finite diquark size, which was very important for the zero density 

case (Fig. 7) to reduce F ln to reasonable values, increases with increasing 

density. That is, the diquark form factor at finite density further suppresses 

the positive contribution of the diquark diagram in Fig.l. The orbital form 

factor Gin shown in the insert of Fig. 11 again demonstrates the enhancement 

of the nuclear magneton, but we have to keep in mind that the backflow effect 

will change the neutron orbital g-factor from the present value to roughly 

| (1 — M m /M N ) < 0, and that effects of the pion cloud around the nucleon 

15 This is in contrast to the quenching of the axial vector coupling constant observed 
in the quark-diquark calculations of Ref.|14| including the axial vector diquark, and in 
hadronic models 
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are known to further enhance the magnitude of the neutron orbital g-factor. 
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Figure 11: Same as Fig. 9 for the form factor F ln (main part) and the neutron 
orbital form factor G^ n defined in Eq. fl4.1|) (insert). 

The medium effects on the neutron form factor i<2 n shown in Fig. 12 are 
qualitatively similar to the proton case of Fig. 10: They decrease the absolute 
value of the "intrinsic" anomalous magnetic moment of the neutron, but when 
combined with the enhancement of the nuclear magneton the neutron spin 
g-factor becomes slightly enhanced in magnitude, as shown in the insert of 
Fig. 12. The magnetic size of the neutron is not changed much, and the total 
medium effects on the neutron spin form factor are very small. 

An interesting common feature of the results shown in Fig. 9-12 is that 
the medium modifications of the orbital and spin form factors (Gl(Q 2 ) and 
Gs{Q 2 )) always decrease with increasing Q 2 . This is consistent with the intu- 
itive expectation that the internal structure of the nucleon at short distances 
is not influenced much by the mean nuclear fields, and again indicates that 
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these form factors reflect the change of the nucleoli current in the medium 
more directly and transparently than the form factors F\, F 2 themselves, or 
any other combination of them. 
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Figure 12: Same as Fig. 9 for the form factor i<2 n (main part) and the neutron 
spin form factor Gs n defined in Eq. (|4.1j) (insert). 



5 Application: The longitudinal nuclear re- 
sponse function 

As an application of the medium effects discussed in the previous section to 
the calculation of nuclear quantities, we consider the longitudinal response 
function for quasielastic electron scattering. Generally, quasielastic processes 
are the ideal places to investigate medium modifications of nucleon form 
factors, because the on-shell kinematics, which is used in the derivation of the 
nucleon form factors, is justified. Recently, experiments on the polarization 
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transfer in proton knock-out reactions have been carried out in the region of 
quasielastic kinematics at large momentum transfers jS], and the results were 
discussed in connection to the predicted medium modifications of nucleon 
form factors [6]. In this work, however, we will consider the simpler case of 
the inclusive quasielastic response function in the region of lower momentum 
transfers. In previous works based on purely hadronic models, it was 
shown that this quantity is quite sensitive to medium modifications of form 
factors, and here we wish to apply our more microscopic quark description 
to the longitudinal response function in nuclear matter. 

The longitudinal response function in isospin symmetric nuclear matter 
is expressed as 



where Hl{u, q) is the 2-point (correlation) function for two external operators 
j%. In the mean field (Hartree) approximation, this is expressed by the first 
diagram of Fig. 13, and if we include the direct terms of the NN interaction 
in the ladder approximation, this gives the familiar RPA series. 



Figure 13: Graphical representation of the longitudinal current-current cor- 
relation function. The solid lines represent nucleons, and the dashed lines the 
interaction between nucleons. The external operator is j%. The first diagram 
gives the Hartree response function, and the others the RPA corrections. 

The relativistic calculation of II ^ follows the lines given in Ref . , 
i.e., the bubble graphs in Fig. 13 consist of the nucleon particle- nucleon hole 




(5.1) 



Tip 
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excitations and the Pauli blocking part of the nucleon particle- antinucleon 
excitations (Z-graphs). The other effects, which are not taken into account 
explicitly in Fig. 13, are summarized in the density dependent electromag- 
netic vertices. Previous calculations based on hadronic models|3*T] incorpo- 
rated the vacuum fluctuations on the level of nucleons (that is, the change of 
the nucleon-antinucleon vacuum polarization graphs in the presence of the 
nuclear mean fields), but it is more appropriate to describe these vacuum 
fluctuations on the level of quarks 16 . Therefore we use our nucleon current 
(|3.1|) for jj, = at the electromagnetic vertices of Fig. 13. 

For the RPA calculation, we need the NN interaction kernel in our ef- 
fective quark theory, which is shown graphically in Fig. 14 and expressed as 

17 



V NN (k) = V a (k) (l)j ■ (1) 2 + V u (k) ( 7a 0i • (Yh (5.2) 

f \ 2 

X 1 F*(k 2 ) (l) x • (1) 2 (5.3) 



-2G V fdMj, N 2 



1 - 2G W U a (k 2 ) + 2G T 5M 2 S V dM 
2C 

Here the first factor in ()5.3|) is the reduced qq t-matrix in the sigma meson 
channel, and the corresponding bubble graph fl a is given in Appendix C. The 
derivation of the effective NN interaction in the NJL model^T], however, has 
shown that in addition to the part Il a , which describes the qq exchange (see 
Fig. 14), there is also a nuclear part which consists of (i) the Z-graph, and 
(ii) a contact term arising from an induced Na 2 N interaction. Concerning 

16 Although this follows naturally from the derivation of the nucleon lagrangian from the 
quark lagrangian in the path integral approach , or from the derivation of the effective 
NN interaction in quark theories following the Landau-Migdal approach 1 1 , it remains to 
be demonstrated explicitly for electromagnetic quantities. 

17 The quantities (1), and (7 M )i in (|5.3(l and l|5.4|l express the Dirac matrices acting 
between the spinors of nucleon i=l,2, i.e., in order to get the NN interaction one has to 
multiply the spinors y/ Mm / 'EmUm of the initial and final states. 
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Figure 14: The NN interaction of Eq. (j5.2j) in the present quark-diquark 
model. Only those terms arising from the interaction of the two specta- 
tor quarks are shown, there are also equivalent graphs from the interaction 
between the quarks inside the diquarks. The term SM 2 , which is included in 
Eq. (J5.3j) . is not shown explicitly in this figure. 

the Z-graph, we note that this is just the Pauli blocking part to the NN 
bubble graph, which is taken into account explicitly in the RPA series of 
Fig. 13 and therefore should not be included in the interaction. (Numerically 
the Z-graph contribution is small compared to n CT because of the reduced 
crN coupling in the medium, see Ref.|ll|.) The self energy correction for the 
sigma meson arising from the induced Na 2 N contact interaction is included 
in (j5.3J) as the density-dependent constant [TT] 

SM ° = i {l^)lTM^) &(VF - M) - (5 ' 5) 

The NNa coupling constant at zero momentum is proportional to dM^/dM, 
and its square appears in (|5.3|) . The NNa vertex form factor F a {k 2 ) is 
normalized to F a (0) = 1. The first factor in ()5.4|) is the reduced qq t- 
matrix in the u;-meson channel, and the corresponding bubble graph fly is 
the same as in the VMD correction to the quark form factors (Fig. 4). The 
NNu coupling constant at zero momentum is proportional to the number 
of quarks in the nucleon, and the vertex form factor F^k 2 ) is normalized as 
i^(0) = 1. 
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Figure 15: The NN interaction for |q| = 0.5 GeV as a function of lo. The 
lower dashed line shows the attractive part V a of Eq. ()5.3|) , the upper dashed 
line the repulsive piece of Eq. (j5.4j) . the solid line shows their naive sum, 
and the dash-dotted line shows the combination V a + V w (l — u 2 /\q\ 2 ). 



By taking matrix elements of IJ5.3J) . ()5.4|) between the nucleon spinors, it is 
easy to see that for k = and for nucleons at the Fermi surface we just get the 
Landau-Migdal interaction derived more generally in Ref . ( |1 lj ) (except for 
the Z-graph contributions as explained above). The form factor F u is equal 
to (Fi p + F ln ), which was calculated in the previous sections. The scalar form 
factor F a should in principle be calculated independently from the Feynman 
diagrams of Fig.l by using the external operator 1. However, because the 
calculations discussed below show that the RPA effects are numerically not 
very important, we will simply assume the same form factor as for the NNu 
coupling (F a = F U ). 

The two parts of the NN interaction, V a of ()5.3|) and V w of (|5.4j) , are 
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shown by the lower and upper dashed lines, respectively, in Fig. 15 for the 
kinematics needed in the calculation of the response function, \q\ = 0.5 GeV 
and < uj < 0.4 GeV. As in the case of relativistic hadronic models [33], 
we see large cancellations between the attractive scalar and repulsive vector 
parts. The solid line shows the naive sum of the two dashed lines, while the 
dash-dotted line shows the combination V a + V w (l — u 2 /\q\ 2 ), which includes 
the effect of the longitudinal space component of the to exchange and is more 
relevant for the longitudinal response function [ST] 18 . The results shown in 
Fig. 15 can be reproduced almost exactly by an approximate form in terms 
of Yukawa potentials, if the coupling constants and meson masses are defined 
at k 2 = 0. This is discussed in Appendix C, where also numerical values are 
given. 

The result for the longitudinal response function of NM (p = 0.16 fm~ 3 ) 
for \q\ = 0.5 GeV is shown in Fig. 16. The dashed line is the Hartree response 
with the free dipole form factors, and the dash-dotted line is obtained by 
adding our calculated medium corrections AF(q 2 ; p) = F(q 2 , p)—F(q 2 , p = 0) 
to the free dipole form factors. (Here F denotes any of the Dirac-Pauli form 
factors.) We see that, even in this region of relatively low momentum trans- 
fers, the medium effects are appreciable. Finally, the solid line shows the 
result obtained by further adding the RPA corrections with the NN inter- 
action derived above and using the density dependent meson-nucleon form 
factors. 

We do not show a comparison to experimental data in Fig. 16 because 

of two reasons: First, our calculation refers to NM and cannot be applied 

18 We have to note, however, that the true interaction in the nuclear medium in the lon- 
gitudinal channel is more repulsive than shown by the dash-dotted line in Fig. 15, because 
of the difference between the Dirac matrices 1 and 7 in Ea. l|5.2|l . On the average it is 
repulsive on the small u> side and becomes attractive on the large uj side, as the RPA result 
of Fig. 16 shows. 
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Figure 16: The longitudinal response function calculated in nuclear matter 
for density p = 0.16 fm~ 3 and \q\ = 500 MeV as a function of u. Dashed 
line: Hartree response with free dipole form factors, dash-dotted line: Hartree 
response including the medium corrections to the nucleon form factors, solid 
line: RPA response including the medium corrections to the nucleon form 
factors. The NN interaction used for the RPA calculation is the one shown 
in Fig. 15, but multiplied by meson-nucleon form factors including medium 
corrections as explained in the main text. 



directly to finite nuclei, although for the case of Ca the results are qual- 
itatively very similar to the NM results, see Ref. j31|. Second, the analysis 
of the experimental data is still controversial, mainly because of the model 
dependence of the Coulomb corrections [S3 IH EH] 19 . 

19 The analysis of Ref.|H], which did not take into account the Coulomb corrections, was 
refined in Ref.0] by using a particular theoretical model for the Coulomb corrections, 
but it has been pointed out recently [201 that the results depend on which theoretical 
prescriptions are used. 
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6 Summary 



In this paper we used a simple quark-scalar diquark picture for the single 
nucleon to describe the electromagnetic form factors of a nucleon bound in 
the nuclear medium. We used the nuclear matter equation of state derived 
within the same effective quark theory to assess the effect of the mean nuclear 
fields on the internal quark structure and the form factors of a bound nucleon, 
taking into account also the meson cloud around the constituent quarks. We 
have shown that this simple model gives reasonable results for the Dirac 
form factors of the free nucleon if the finite size of the diquark is taken into 
account. This is particularly important for the neutron in order to obtain a 
small Fi, consistent with observations. Concerning the Pauli form factors, in 
particular the anomalous magnetic moments and the behavior for large Q 2 , 
one would need to further add the effect of the axial vector diquark and its 
finite size to achieve a reasonable description. 

The medium modifications of the form factors associated with the orbital 
current are significant in the region of low and intermediate Q 2 , and partially 
associated with an increase of the electric size in the medium. The form 
factors associated with the spin current are enhanced because of the reduced 
nucleon mass, but due to a simultaneous decrease of the intrinsic anomalous 
magnetic moment, the total changes are not very large; in particular for 
the neutron they are almost zero. For both kinds of form factors - orbital 
and spin -, the medium modifications decrease with increasing Q 2 . This is 
consistent with the intuitive expectation that the mean fields, which reflect 
the long range nuclear correlations, should not influence the structure of the 
nucleon at short distances. 

As an application, we considered the longitudinal response function of 
nuclear matter for inclusive quasielastic electron scattering. Even in this re- 
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gion of relatively low momentum transfer, the medium effects on the nucleon 
form factors give rise to an appreciable quenching of the response function. In 
order to take into account the RPA-type correlations between the nucleons, 
we derived the NN interaction in the present quark model for the isoscalar 
channel, which is most important for the longitudinal response function. In 
the limit of zero momentum transfer (Landau limit), this interaction agrees 
with the more general Landau-Migdal effective interaction, which was de- 
rived in earlier works for this particular quark theory, and has many features 
which are similar to relativistic meson-nucleon theories. In particular, it is on 
the average repulsive in the region below the quasielastic peak, and becomes 
attractive for higher energy transfers. 

We finally would like to remark that the language of effective quark the- 
ories is very appropriate to address the problem of medium modifications of 
nucleon properties. This work represents one further step toward the goal 
of describing relativistic nuclear systems by taking into account the quark 
substructure of the constituent hadrons. 
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Appendices 



A Expressions for the nucleon form factors 

Here we give the explicit expressions for the contributions of the 3 terms in 
Eq. (|3.11|) - (j3.13|) to the nucleon Dirac-Pauli form factors defined in (|3.14|) . 
They can be derived by invariant integration in the usual way, i.e., by (i) 
introducing Feynman parameters, (ii) performing a shift so that the denom- 
inator depends only on k 2 , (iii) expressing the result in the form (j3.14|) by 
using the Dirac equation for the nucleon spinors, (iv) performing a Wick 
rotation according to 

d4 W 2 )-fS/H) 



(2tt) 4 '' v ' Jo 16vr 2 ' 
where t = k 2 + k 2 is the square of the Euclidean length, and finally (v) 
introducing a Lorentz invariant regularization scheme, in our case the proper 
time regularization (j2.19j) . Below we give the expressions which are obtained 
after performing the steps (i)-(iv), using Q 2 = —q 2 > as the variable. 

In the expressions given below there enter the diquark and nucleon wave 
function normalization factors, which are defined by (j2.13J) and (j2.14j) in 
terms of the renormalized bubble graphs (j2.15|) and (j2.16|) . for which we 
have the expressions 



oo 



, tdt ( i p 2 r 1 , i 

ILV) = -24 / — - — — + 1- dx 



o 16tt 2 \t + M 2 2 Jo (t + M 2 -p 2 x(l - x)f ) 

(A.l) 

r°° tdt f l p'x + M 

Rn{P) = -g, / —r^ / dx 



If we denote 



o 16Wo (t + M 2 {l-x) + M 2 s x-p 2 x{l-x)f 

(A.2) 



D!(Q 2 ,x) = t + M 2 + Q 2 x(l - x) , (A.3) 
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the contributions of the contact term (J3.11|) to the form factors are as follows: 

Ff N \Q 2 ) = 



F 2N (Q 



If we denote 
D 2 (Q 2 ,x,y) 

m(Q 2 ,x) 

N 2 {Q\x) 



AG S ~ r°° tdt yi 

— Z N / T7j 2 / dx 

g s Jo Ibir z Jo 



Fiq(Q 2 



1 



^!(Q 2 ,x) t + M 2 2 D^Q^x 



+ F 2 q(Q 



Q 2 (l+x)(l-x 



M - 



2 A(Q 2 ,x) 



(A.4) 



^z N M N Mr^Lrax 

g s Jo Ibn 2 Jo 



Fiq(Q 2 



Di(Q 2 ,x) 



(-1 _ ^Mjv\2 _ Q^x 

+ F2Q{Q } D 1 (Q 2 ,x) 2 



(A.5) 



O 2 

t + M 2 (l - x) + M 2 x - M 2 N x(l -x) + -j- {x 2 - y 2 ) 



(M N + M) 
M N 



2 -Ml-*-- 2M 2 N x{l -x)- 2M N Mx 



1 + 



M 



AM 2 



(x 2 - y 2 ) 



(A.6) 



the contributions of the quark term (j3.12j) to the form factors follows: 



IN 



Jo lbir z Jo 

W 5! r,.>±f(i- 



* N x {Q 2 ,x) 



D 2 (Q 2 : x,y)\ 



dy 



7 (Q) 

IN 



'N 



2 J- x ~" D 2 (Q 2 ,x,y)* 

tdt rl 



(A.7) 



167T 2 Jo 

F 1Q (Q 2 )2M N x 



dx / dy 



, n . r M N (l-x) + M , „ (ra . % , „ N 2 (Q 2 ,x) 



D 2 (Q 2 ,x,yf 



+ F 2 q(Q ) M N M 



D 2 (Q 2 ,x,y)3 
(A 
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If we denote 

2 

D 3 (Q 2 , x,y) = t + M 2 (l-x) + M 2 x - M 2 N x{l - x) + -^-(x 2 - y 2 ) , 

(A.9) 

the contributions of the diquark term (j3.13j) to the form factors are as follows: 

'°° tdt /"i /•* 2M N (1 - x)(M N x + M) + | 



/o 167r 2 jo J-x D 3 (Q 2 ,x,y) 3 

(A.10) 

(A.ll) 

The diquark form factor Fp is calculated from the expression ()3.16|) as 

poo fAf pi ( . . | ^a; x 

where we defined 

£> 4 (Q 2 , x,y) —t + M 2 — M 2 x(l - x ) + ^-(x 2 - y 2 ) , (A.13) 
and F^q , F^q are the isoscalar parts of the quark form factors. 

B Expressions for the quark form factors 

Here we give the expressions for the quark electromagnetic vertex 1)3.25)1 in 
terms of the quark form factors defined in Eq. (J3.8j) . In the expressions given 
below there enter the pion and quark wave function normalizations g n and 
Zq, which are defined in terms of the renormalized self energies 11^ and Eq 
by Eqs. (|3.24j) and ()3.26|) . The expression for iX = II S has been given in 
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Eq.flHH), and 

a ( , f°° tdt /•! ^x-M 



o 16tt 2 7o (t + M 2 (l-x) + M 2 x-p 2 x(l-x);r 

(B.l) 

As explained in the main text, the isoscalar (or isovector) parts of the quark 
vertices given below should eventually be further multiplied by the VMD 
form factors ()3.32|) . where the expression for fly is 

The quark diagram (second term in (|3.25Jl ) gives the following contribu- 
tions to the quark form factors: 

(B.3) 

r,(Q) 1 /, \ ~ /*°° Mt r 1 , /■* 2M 2 x 2 



where 



D 5 (Q 2 , s, i/) = t + M 2 (l - x) + MV + ^(x 2 - y 2 ) . (B.5) 

The pion diagram (third term in (|3.25jl ) gives the following contributions 
to the quark form factors: 



F$ = -2r 3 g w F,(q 2 ) I / dx / dy 



°° tdt /-i r x , 2M 2 (1 -x) 2 - I 



/o 167r 2 Jo j-x D 6 (Q 2 ,x,y) 3 

(B.6) 

4? - 4 -M 2 ^^ 2 )/;^f dx£d,^^, 

(B.7) 



where 



O 2 

D 6 (Q 2 ,x,y) = t + M 2 (l-x) 2 + M 2 x+^-(x 2 -?/ 2 ). (B.8) 
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The pion form factor is calculated from (j3.23j) and (j3.28j) as 

poo pi | px 



x 



D 7 (Q 2 ,x, y y 



where 



D 7 (Q 2 , X ,y)= t + M 2 - M 2 x(l - X ) + ^(x 2 - y 2 ) 



C The NN interaction 



(B.9) 



(B.10) 



Here we provide some details on the NN interaction Eqs. fj5.2jl . 

The bubble graph in the sigma channel, which appears in (|5.3j) . has the 
form 



6/ 
12 



d 4 g 
00 tdt 



Ti D (S F (q)S F (k + q)) 
2 



167T 2 



t + M 2 



+ (k 2 - AM 2 ) J dx 



1 



(t + M 2 -k 2 x{l-x)f 

(C.l) 



The expression for Uy has been given in (jB.2j) . 

In the calculations of Sect. 5, the NN interaction (j5.2j) is used without 
further approximations. Here we wish to give an approximate form in terms 
of Yukawa potentials: Expanding Il a and fly around k 2 = 0, we obtain 



Vnn(&) 



9<tN 



k 2 -Ml 



'l)-(l) 



9u>N 



k 2 - M 2 



(7,) • (v) 



(C.2) 



where 



'dM 



9uN 
9luN 



N 



g a = (2.34) 2 x 21.75 = 119.4 



V dM f 
9g u = 9 x 18.43 = 165.9 



(C.3) 
(C.4) 
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, and the numerical 



where the quark-meson couplings are defined by 
g a = [(dfljdk*) k2= y and g u = [- (dfl v /dk>) ^ 
values given above are obtained from our nuclear matter EOS for p = 0.16 
fm~ 3 . The meson masses defined at zero momentum are 



Ml 



Mr. 



9a 

9u 



1 



2G n 



2a 



U a (0) + 5Ml) = (0.81 GeV) 1 
;i.l3GeV) 2 . 



(C.5) 
(C.6) 



Note that these masses are different from the pole positions. 
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